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ABSTRACT 

The force on an electron in an ionized medium may be written (ff+^iraP per 
unit charge, where & is the electric force, P the polarization of the medium, and a 
a constant which is zero on the Sellmeyer theory and may have a value of 1/3 on 
the Lorentz theory. A well-defined distinction exists between maximum usable 
frequencies for ionospheric radio transmission calculated on the basis of the two 
theories. This paper describes a type of experiment for measuring maximum 
usable frequencies and comparing them with calculated values. The results 
indicate that a is probably zero, at least in these experiments. There is also 
described a stQp-by-step method of solving the virtual-height integral equation, 
necessary in much theoretical work concerning the ionosphere. 
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I. INTRODUCTION 

The question whether the force on an electron in an ionized medium 
should be taken as merely the applied electric force, fi } , or the so-called 
Lorentz force, $ +47raP [l], 1 where P is the polarization of the medium 
and a is a constant which is usually taken as 1/3, is of fundamental 
importance in the theory of the ionosphere. The two theories are 
known, respectively, as the Sellmeyer theory and the Lorentz theory, 
and the difference between them represents a difference of 50 percent 
in the equivalent electron density in the ionosphere. 

The argument had apparently been settled in favor of the Sellmeyer 
theory in a comprehensive theoretical paper by Darwin [2], when the 
whole question was reopened by a radio experiment of Booker and 
Berkner [3], which they were able to interpret only in terms of the 
Lorentz theory. In Australia, Martyn and Munro [4] explained the 
same kind of radio experiment in terms of the Sellmeyer theory, but 
Appleton, Farmer, and Katcliffe [5] in England attacked this inter- 
pretation. 

The radio experiments referred to were all performed at vertical 
incidence, in the neighborhood of the gyrofrequency, or precessional 
frequency of electrons in the earth's magnetic field (about 1.4 Mc/s). 
The purpose of this paper is to present the results of a different kind 

» Figures in brackets indicate the literature references at the end of this paper. 
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of radio experiment — one involving ionospheric transmission of radio 
waves over a distance at oblique incidence, at a high frequency, and 
during the daylight hours. 
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Figure 1. — Oblique-incidence relations for Sellmeyer and Lorentz theories. 

z/=height of equivalent triangular path, a=l/3 
z,=height of equivalent triangular path, a=0 (equal to virtual height at equivalent vertical-incidence 

frequency) 
/=equivalent vertical-incidence frequency, a =1/3 
/' cos 0i=equivalent vertical-incidence frequency, a=0 

0i= angle of incidence of waves upon the ionized layers. 

II. THEORY, NEGLECTING THE EARTH'S MAGNETIC 

FIELD 

There exists a well-defined distinction in the behavior of oblique- 
incidence radio transmission according to the two theories. This 
distinction is that the maximum usable frequency for radio sky-wave 
transmission over a given distance is greater for the Lorentz theory 
than for the Sellmeyer theory. 

In two previous papers [6, 7] the author outlined a method of calcu- 
lating maximum usable frequencies and heights of reflection for 
oblique-incidence radio transmission directly from the curve of 
frequency versus virtual height obtained at vertical incidence. This 
curve will hereafter be called the (z v J) curve, where z v is the virtual 
height and f is the vertical-incidence frequency. The virtual height 
is obtained by measuring the time, M f taken by a pulse to travel up 
to the ionosphere and back. It is given by 
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where c= velocity of light, 

/= frequency of waves, 
ju= refractive index of the ionosphere, 
z = actual level of reflection, and 
z= measured vertically upward. 
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The method of calculating maximum usable frequencies was based 
on the simple ray theory of a stratified isotropic ionosphere, whose 
refractive index, /x r , for waves of frequency/' was given by 



: y l ~p 



a) 



where/' is the actual wave frequency and j 2 —Ne 2 /irm. } N being the 
ionization density expressed as a function of height. 

The calculation of maximum usable frequencies consists essentially 
in solving graphically the two equations z v =z v (j) 1 the (z v J) curve for 
vertical incidence, and s«=2//,jf / , D), the transmission equation for 
oblique incidence, obtained below. Here D is the distance of trans- 
mission,/' the wave frequency, and/ the equivalent vertical-incidence 
frequency, or frequency of the wave reflected, at vertical incidence, 
at the same level as is the wave frequency, /', over the given dis- 
tance, D. 

The transmission equation is found by combining the expression for 
the equivalent triangular path, obtained geometrically, with the 
expression for the equivalent vertical-incidence frequency, /, and the 
relation between z v and z/, the height of the equivalent triangular 
path. 

The intersection of the graphs of the two equations gives both the 
equivalent vertical-incidence frequency, /, and the virtual height of 
reflection corresponding to transmission of the given wave frequency, 
/', over the given distance, D. The greatest value of/' for which a 
real solution exists is the maximum usable frequency over the given 
distance. 

When the Lorentz polarization term a is included, the refractive 
index, //, is given [8] by 



"'"V 1 -/^- (2) 

The virtual height at the vertical-incidence frequency, /, is 

Jo M 

where z is the level of reflection, and \x is n' with/ substituted for /'. 
The equivalent vertical-incidence frequency, /, is given by 



^Vl-acos 2 ^ 



/=/' cos <j>i\ i J-T* (4) 

where <j>i is the angle of incidence of the actual wave upon the iono- 
sphere. 

The height, z v ', of the equivalent triangular path may be calculated 
as follows. Let <t> be the angle between the ray path and the vertical 
at the height z. Then 
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dz 



J*o ^ z C z o = 
- = COS 0i I / o , fo- 
il' COS I - /COS 2 0i — 



f' 2 +af 2 



This becomes, on substituting/ for/', 



-XV 



(1 — a) cos 2 0^ 2 
1-acos 2 */ ^ (5) 



f-(l-a)f 2 
To obtain the transmission equation ,we get geometrically the value of 

We then insert in eq 6 the value of t obtained from eq 4 in terms of/ 
and/', and multiply 2/ by the ratio z v : z v ' determined by comparing 
eq 3 and 5. 

This analysis has been made for a flat ionosphere, but may be 
extended to a curved ionosphere in the manner indicated in reference [7]. 

For a=l/3, the ratio f j1 is greater than for a=0, and so the wave 
frequency,/', corresponding to a given/ and z v is in general greater. 
If z v were equal to z/ f therefore, the calculated maximum usable 
frequency would be much greater for a=l/3 than for a=0, the dif- 
ference increasing with the angle of incidence, $ la We must, however, 
examine the relation between z v and 2/. Now z is a function of / 2 , 
so if we put x=f 2 /f 2 we may rewrite eq 3 and 5 : 
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Equation 7 gives z 9 as a function of/ 2 and eq 8 gives 2/ as a function 
of/ 2 and cos 0i, so that z v /z v ' may be expressed as a function of/ 2 
and cos 0i. 

The ratio z v [z/ is, for a=l/3 ) in general less than unity, that is, less 
than the value it had for a=0. This tends to decrease the calculated 
maximum usable frequencies. To what extent the decrease due to this 
effect compensates for the increase due to the greater ratio of/'// must 
be determined in individual cases. 

Let us first consider the case of a linear gradient of ionization 
density, where z=kf 2 . Here dz/dx=kf 2 and the determination of 
z v and 2/ is easy. We obtain 

^=K^)[- i+ v^ tan_l VS] <» 

g /= i"l + , * m tan-' J J cos2 *' -1. (10) 

L V^ cos 2 0i(l— a cos 2 00 V 1— a cos 2 0! J 
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Figure 1 shows the ratio of z t '/z 9 for a=l/3, which is the same as the 
ratio of z/ for a=l/3 to zj for a=0, plotted against cos (fo. It shows 
also the ratio (/'cos 0i//) 2 for a=l/3, which is the same as the square of 
the ratio of/ for a=0 to /for a= 1/3, plotted against cos 0j. The two 
ratios are equal at about cos 1 =O.8. This means that, for the 
distance corresponding to this value of cos b the wave frequency 
corresponding to reflection at a point on the linear part of a (z v J 2 ) 
curve is the same for a—0 as for a=l/3. At shorter distances the 
wave frequency is less, and at greater distances it is greater, for a— 1/3 
than for a—0. 

If the (z v J 2 ) curve curls up, as happens near a critical frequency, 
an increase in the height of the equivalent triangular path has rela- 
tively less effect than in the case of a linear (z„J 2 ) curve. This means 
a relative increase in the wave frequencies calculated for points on a 
curve near a critical frequency. Thus it may be expected that the 
maximum usable frequencies calculated for an actual distribution of 
ionization density will be, for a=l/3, equal to or less than those cal- 
culated for a=0 only for very short distances. Elsewhere they will 
be greater, and may become as much as 20 percent or so greater, at 
the greater distances, than in the case where a—Q. The calculation 
must in general be made for each individual distribution of ionization 
density. 

The calculation can be made for a general distribution of ioniza- 
tion by assuming the (z v J 2 ) curve to be made up of a number, m, of 
linear intervals. In the nth. such interval the ionization density, 
represented by/ 2 , is assumed to vary between /„ 2 and/ n+1 2 , while the 
height varies from to z n . In this interval, then 

*\JQ J / 2 / 2 z n* 

Jn Jn—1 

Then, if x=j */f and/ m 2 =/, 

dz_ fz„ 



dx / n+1 2 -/„ 2 
and is constant over the interval, so that eq 7 becomes 
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Each of the integrals in the summation represents a part of the 
complete integral from to 1/(1— a) and may be readily evaluated 
for a given a and cos 4>i- Both eq 11 and 12 are linear in z n and may 
easily be evaluated. If we know only the (z v J) curve, eq 11 may be 
turned inside out and the known values of z v used to form n equations 
in the n unknown values of z n , which can then be substituted in eq 12 
to obtain zJ. 



110 Journal oj Research oj the National Bureau oj Standards [Vol. $6 



Calculations of the maximum usable frequencies for E- and F 2 -layer 
transmission were made, using this method, for a=0 and a=l/3, for 
typical distributions of ionization density. The ratios of the maximum 
usable frequency for a=l/3 to that for a=0 are plotted, for various 
values of distance D, in figure 2. The solid graph is for i^>-layer and 
the dotted for Z?-layer transmission. 

For i^-layer transmission, it is seen that the maximum usable 
frequencies are not much greater for a=l/3 than for a=0 until the 
distance of transmission exceeds about 500 kilometers. For greater 
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Figure 2. — Ratios of maximum usable frequencies for a- 
E- and F 2 -layer transmission. 



1/3 to those for a=0 f for 



distances the maximum usable frequencies for a— 1/3 are higher than 
those for a=0, becoming about 17 percent higher at the greater 
distances. 

For J^-layer transmission the ratio of maximum usable frequency 
for a=l/S to that for a=0 rises much more rapidly with distance than 
in the case of 2^-layer transmission, and attains a somewhat greater 
value. 



III. EFFECT OF THE EARTH'S MAGNETIC FIELD 

The calculations so far have not included the effect of the earth's 
magnetic field. Since this field may affect considerably some high- 
frequency transmissions, it is necessary to consider its effect. Accord- 
ingly calculations were made for the transmission path and frequency, 
for which experimental results are given below. 

The transmission path was practically perpendicular to the magnetic 
meridian, and was about 650 km long. The earth's field was assumed 
to have an intensity of 0.5 gauss, and a dip of 70°. For a frequency of 
6.06 Mc/s, and a vertical-incidence virtual height of 130 km at the 
equivalent vertical-incidence frequency, table 1 gives the calculated 
oblique-incidence data. These were calculated on the basis of an 
average distribution of ionization density in the E region, using the 
magneto-ionic equations as published by Booker [91. 
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Table 1. — Calculated oblique-incidence data 
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2.31 
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The table gives, for a=0 and 1/3 and for 11=0 and 0.5 gauss, the 
ionization density at the level of reflection, the angle of incidence of the 
waves which traverse the given path, the equivalent vertical-incidence 
frequency, and the vertical-incidence critical frequency which cor- 
responds to a maximum usable frequency of 6.06 Mc/s over the given 
path, via Z?-layer transmission. It may be seen that the inclusion of 
the earth's field where a=0 makes a difference of but 2.0 percent, and 
where a=l/3 a difference of but 1.0 percent. The inclusion of the 
polarization term, however, makes a difference of 15.1 percent where 
H=0 and 14.2 percent where 11=0.5 gauss. 

Thus it should prove easy to decide between the Lorentz and 
Sellmeyer theories by measuring the vertical-incidence i?-layer 
critical frequency at the times when the oblique-incidence E'-layer 
transmission begins and fails. The experimental error should be con- 
siderably less than the difference between the two theories, either for 
calculations based on H=0 or for closer calculations based on 11=0.5 
gauss. The calculations for east-west transmission may be made for 
H=0 with an error much less than the probable experimental error. 

A recent analysis has been published by Ratcliffe [10] for a para- 
bolic distribution of ionization density in the ionosphere. The curves 
and discussion he gave for oblique-incidence radio transmission indi- 
cated little difference between the results for a=0 and a=l/3, es- 
pecially for long-distance transmission via the F layer. Indeed, he 
stated that the difference would be greatest for a transmission distance 
of about 500 km. 

The curves of figure 2 make this statement appear in error. It 
seems likely that Ratcliffe was misled by making calculations only for 
ranges below about 700 km, as published in his paper, for it is obvious 
that the greater separation between the two cases comes for consid- 
erably greater distances. It is possible, moreover, that the assump- 
tion of a parabolic distribution might also tend to lessen the difference 
at the shorter distances. 

The results of the step-by-step analysis of typical E- and i^-layer 
curves indicate the nature of the difference in maximum usable fre- 
quency according to the two theories. 

IV. EXPERIMENTAL OBSERVATIONS 

There is thus a clear distinction between the Lorentz theory 
(a=l/3) and the Sellmeyer theory (&=0) in the case of radio waves 
incident obliquely upon the ionosphere. It should therefore be pos- 
sible to decide experimentally between the two theories. 
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There are available considerable published data on oblique-incidence 
radio transmission [11, 12, 13, 14, 15]. Experiments involving the 
E layer, which, as shown by figure 2, should provide conclusive 
results for relatively short transmission distances, are for the most 
part unreliable, usually because of the presence of sporadic-Z? reflec- 
tions above the true 25-layer critical frequency. Such an experiment 
was attempted by Farmer and Ratcliffe [11], with somewhat incon- 
clusive results. Experiments involving the F 2 layer have for the 
most part been conducted over distances of less than 600 kilometers. 
They are subject to the difficulty that ionosphere conditions at the 
point of reflection are much more uncertain for the F layer than for 
the E layer. Also, the difference between the two theories would 
produce a discrepancy of but a few percent, which is of the same 
order of magnitude as the uncertainty of ionosphere conditions at 
the point of reflection. Much of the published data on radio trans- 
mission, therefore, is of little use in deciding whether a=0 or 1/3. 

Since September 1935, the National Bureau of Standards has made 
continuous automatic field-intensity records of numerous radio sta- 
tions. Results are available for the whole period from 1935 to 1940 
for high-frequency international broadcasting station W8XAL (later 
WLWO), 6.06 Mc/s, at Mason, Ohio, 645 kilometers distant. Some 
results of these recordings were presented in a paper before the joint 
meeting of the International Scientific Radio Union and the Institute 
of Radio Engineers at Washington in April 1936, and other results 
have been used by investigators at the Bureau in studies of the 
ionosphere [16, 17]. 

The transmissions from W8XAL to Washington were propagated 
normally by way of the E layer during most of the daylight hours, 
and by way of the F layer at other times. The transitions from F- to 
j^-layer transmission in the morning, and vice versa in the afternoon, 
were often well marked, since the vertical radiation pattern of the 
transmitting antenna was such as to provide substantially greater 
radiation for single-hop E- than for i^-layer transmission. The 
transitions have been identified on a great many records by the 
following facts. (1) The regularity of occurrence of the change in 
intensity and character of fading on the records corresponds to the 
regularity of the diurnal and seasonal variations of the i?-layer 
critical frequency. (2) There was a complete absence of any regular 
sudden change in local ionosphere characteristics (critical frequency, 
virtual height, or absorption), except for the jB-layer critical fre- 
quency, at the time of the change on the record. (3) The beginning 
of propagation was observed when no other ionosphere layer but the 
E layer had sufficient ionization density to support transmission (as 
during the mornings of ionosphere-storm days). (4) The change on 
the record nearly always occurred about 15 minutes after the local 
Z?-layer critical frequency rose or fell through about 2.45 Mc/s. This 
time difference is the difference between the local times at Washington 
and at the midpoint of the path from W8XAL to Washington. 

Figure 3 shows four examples of the change of layer on the W8XAL 
records. The record of May 12, 1938, was made when a severe 
ionosphere storm had considerably depressed the critical frequencies 
of the F x and F 2 layers. At the time when good transmission began, 
none of the layers except the E layer had even nearly enough ioniza- 
tion to support transmission. Recordsof February 21 and22, 1938, were 
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Figure 3. — Transitions between E- and F 2 'layer transmission from Mason, Ohio, to Washington, D. C. 

Station W8XAL, 6.06 Mc/s, 650 km distant. 
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made during a period when the night and morning i^-layer transmis- 
sions from W8XAL were greatly absorbed. The transitions here are 
seen to be quite sharply defined and symmetrical about local noon at 
the midpoint of the transmission path. The record of March 13, 1937, 
is typical of most of the normal W8XAL records on which the transi- 
tion can be seen. The transitions are not quite so sharp, but are 
easily identifiable, and-jthe times may be measured within a few min- 
utes. In the files of the National Bureau of Standards are 11 records 
similar to that of May 12, 1938, 13 similar to those of February 21 
and 22, 1938, and several hundred showing the transition like that of 
March 13, 1937. 

For the average regular distribution of ionization in the E layer, 
a maximum usable frequency of 6.06 Mc/s over a distance of 645 km 
corresponds to a vertical incidence/* of about 2.40 Mc/s for a=0, 
and 2.10 Mc/s for a=l/3. The midpoint of the great-circle path 
between jW8XAL and Washington is almost due west of Washington 
and the distance is short enough so that the/* may be considered a 
function of local time only. Therefore the W8XAL transmissions are 
well adapted for deciding the value of a. 

Since the E layer at the time of these observations is known to be 
chiefly electronic in character, it would seem that the presence of 
heavy ions cannot be invoked to explain a value of a different from 
1/3 or 0. Also, it should be noted that it takes about 32 minutes for 
ihef B ° to rise in the morning, or fall in the evening, between 2.06 and 
2.40 Mc/s. Also, the E x reflections were so absorbed at these times, 
on the records analyzed, that there can be no question of confusing 
E°- with E'Mayer transmission. 

Table 2 shows the result of comparison of 13 records, in each of 
which at least one transition was as well-defined as in the record of 
February 21 or 22 of figure 3. The values of a corresponding to the 
observed ratio of maximum usable frequency to vertical-incidence 
critical frequency are tabulated. In calculating these values of a, 
an average distribution of jE-layer ionization was used instead of the 
actual distribution for each day, which varied somewhat from day to 
day. This may account for some of the variations in a. 

The average value of a deduced from the 19 transitions listed in 
table 2 is —0.02, or —1/50, as against the theoretical value of 1/3. 
This value of a is not believed to be significantly different from zero. 
A number of ionosphere storm days were studied, when no layer 
other than the E layer was sufficiently ionized, in the morning, to 
produce transmission. The transitions on many of these days were 
obscured by sporadic-2? reflections. The mean value of a calculated 
on 7 days when the transition could be easily identified was —0.005, 
again not significantly different from zero. 

Because of the great number of the records showing well-marked 
transitions which were not quite so clear as those just mentioned, the 
transitions on these records were studied statistically, rather than 
individually. The monthly averages of the morning and evening 
transition times were obtained, and the monthly average/* at this 
time was determined from the vertical-incidence ionosphere data. 
Table 3 shows the average values of f E ° at the midpoint of the trans- 
mission path, at the times when JS'-layer transmission began and 
failed, and the corresponding values of a. The mean value of a, 
deduced from 435 transitions studied in this group, was —0.06, 
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again not significantly different from zero. The negative value of a 
probably means that the height of the E layer was slightly greater 
than that assumed as an average. 



Table 2.- 



- Values of "a" calculated from value of local f E at time of change of layer for 
W8XAL, 6.06 Mc/s, 650 km distant 



Date 

(February 

1938) 


Local /b in Mc/s 


Calculated value of a 


A.M. 


P.M. 


A. M.- 


P.M. 


6 
14 
15 
17 
18 
21 
22 
23 
24 
25 
26 
27 
28 


~~2.~33~~ 
2.48 

~~2.~35~" 
2.44 
2.46 
2.48 
2.48 
2.54 
2.39 
2.36 
2.44 


2.47 
2.43 

~"2.~46" 

~~2.~45~~ 
2.48 

~~2.~37"~ 
2.29 
2.35 




-0.08 
-.03 


-4-0.10 
-.10 


-.07 


+.08 
-.05 
-.07 
-.10 
-.10 
-.19 
+.02 
+.06 
-.05 


-.06 
-.10 




+.07 
+.15 
+.08 








Table 3. — Mean values of "a" calculated from values of local fs at time of change of 
layer for W8XAL, 6.06 Mc/s, 650 km distant 



Month 


Morning transition 


Evening transition 


Time 


fs° 


a 


Time 


fit° 


a 


1936 
May 








1823 
1838 
1830 
1737 
1651 
1610 
1559 

1709 

1739 
1759 
1818 
1849 
1852 


2.50 
2.45 
2.42 
2.40 
2.42 
2.47 
2.53 

2.42 
2.40 
2.50 
2.48 
2.40 
2.44 


-0.13 
-.06 
-.02 
+.01 
-.02 
-.08 
-.17 

-.02 
+.01 
-.13 
-.10 
+.01 
-.05 


June 








July 








September 


0703 
0738 
0810 
0831 

0805 
0726 
0637 
0604 


2.44 
2.52 
2.40 
2.48 

2.45 
2.47 
2.40 
2.42 


-0.05 
-.16 
+.01 
-.10 

-.06 
-.08 
+.01 
-.02 


October 


November. 


December.. 


1937 
February 


March 


April. 


May 


June. . 


July... 








Mean a 












-0.06 






-0.06 













V. CONCLUSIONS 

For the cases considered in this paper, the data indicate that the 
Lorentz polarization term should not be included in the theory of 
oblique-incidence radio transmission by way of the ionosphere. Fur- 
ther experiments in different latitudes and with careful determina- 
tions of oblique and vertical-incidence penetration frequencies should 
be performed before this conclusion can be positively confirmed. 

It is suggested that experiments of this type, involving wide ranges 
of frequencies and reflection heights, are less conducive to misleading 
results and ambiguous interpretations than are vertical-incidence 
experiments in the neighborhood of one frequency, which involve 
somewhat uncertain factors, such as the magnitude and direction of 
the earth's magnetic field at definite heights, and assumptions as to 
the physical constitution of the ionosphere. 
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